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Abstract 



We explore how nuclear modifications to the nucleon parton distributions 
affect production of high transverse momentum hadrons in deuteron-nucleus 
collisions. We calculate the charged hadron spectra to leading order using 
standard fragmentation functions and shadowing parameterizations. We ob- 
tain the d+Au to pp ratio both in minimum bias collisions and as a func- 
tion of centrality. The minimum bias results agree reasonably well with the 
BRAHMS data while the calculated centrality dependence underestimates the 
data and is a stronger function of px than the data indicate. 



1 



Typeset using REVTgX 



One of the most intriguing results from the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory has been the suppression of hadrons with high transverse 



momentum, pr, in central Au+Au collisions at center-of-mass energies, y/Srm, of 130 and 
200 GeV. The AA suppression factor, 

c / \ da AA /dp T 

Raa{Pt) = , n rp v , 7^— » I 1 ) 

{a^ N T AA )da pp /dp T 

compares the AA and pp pr distributions of hadrons, normalized by the number of binary 
collisions, (ff^T^), the product of the nucleon-nucleon inelastic cross section, <tnn, and 
the nuclear overlap function T AA . Saturation effects in the initial nuclear wave function and 
final-state parton energy loss were both proposed as explanations of the large suppression 
seen in Au+Au collisions by PHENIX [1,2], STAR [3,4] and BRAHMS [5]. To determine 



whether the suppression is an initial or final-state effect, d+Au collisions at ^/Snn = 200 
GeV were recently studied at RHIC. The data [6-8] show that, at midrapidity (rj ks 0), the 
d+Au suppression factor, 

p ( n \ da dA u/dp T fo s 

{cr^ N T dAu )dap P /dp T 

is much closer to unity. These results suggest that the strong suppression in Au+Au collisions 
is a final-state effect, implying that, at least at central rapidities, saturation effects are small. 
However, at higher rapidities where the nuclear parton momentum fraction, x 2 , is smaller, 
such effects might still be important. Since x is not very small at RHIC, it is necessary to 
check if other, more conventional, models of nuclear shadowing may also explain the data. 

The BRAHMS collaboration has measured -RdAu at several values of pseudorapidity, 77, 
and observed increasing suppression as rj increases from \q\ < 0.2 to rj = 3.2 [8]. The 
BRAHMS measurements are in four 77 bins: < 0.2; 0.8 < 77 < 1.2 (77 = 1); 1.9 < 77 < 2.35 
(77 = 2) and 2.9 < r/ < 3.5 (r/ = 3.2), corresponding to center-of-mass scattering angles, 9 cm , 
of 101.4° > 9 cm > 78.6°, 48.4° > 9 cm > 33.5°, 17.01° > 9 cm > 10.9° and 6.3° > 9 cm > 3.5° 
respectively. These data have also been divided into three centrality bins: (0 — 20)%, 
(30 — 50)% and (60 — 80)% of the geometric cross section. Using a Glauber calculation of 



the nuclear overlap with a Woods-Saxon density distribution for the gold nucleus and the 
Hulthen wavefunction to calculate the deuteron density, we find that these bins correspond to 
the impact parameter, b, ranges of < b < 3.81 fm, 4.66 < b < 6.01 fm and 6.59 < b < 7.74 
fm respectively. Our calculated values of (cJvAr^dAu) in these impact parameter bins are 
in relatively good agreement with those determined by BRAHMS. The results in the two 
lowest rj bins are reported for {h + + h~)/2 while the rj = 2.2 and 3.2 bins are reported for h~ 
only where h + and h~ stand for the positively and negatively charged hadrons respectively. 

In this paper, we calculate RdAuipr) in the BRAHMS i] bins using two parameterizations 
of nuclear shadowing. We also calculate the central-to-peripheral ratios, i?cp(pr), with 
two parameterizations of the spatial dependence of shadowing. The calculated ratios are 
compared to the BRAHMS data [8]. To better illustrate the effects of shadowing alone, we 
do not include the Cronin effect, pr broadening [9,10], in our calculations. 

We make a leading order (LO) calculation of minijet production to obtain the yield of 
high-p r partons [11]. The p T distribution of Ref. [12] is modified to include the nuclear 
parton distribution functions, 

" X = 2pT y~] [ -pp— [ dxi ( dx 2 [ d 2 s [ dz [ dz' 

„• ~~f^ „ JOmin sin 9 cm J J J J J 



j^D h / k (z c ,Q 2 ) diJij^k 



x F i/d ( Xl , Q\ s, z)F j/Au (x 2 , Q\ \b - s\,z') n/K " ^ (3) 

Zq ux 

where x 1 and x 2 are the parton momentum fractions in the deuterium and gold nuclei 
respectively, Q is the momentum scale of the hard interaction and z c is the fraction of the 
parton momentum transferred to the final-state hadron. The integrals over center-of-mass 
scattering angle, ^min < 9 cm < # max , correspond to the BRAHMS angular regions, given 
previously. The 2-^2 minijet cross sections, d&ij^k/di, are given in Ref. [13]. Even though 
the next-to-leading order corrections may affect the shape of the pt distributions, the higher- 
order corrections should largely cancel out in -RdAu, as is the case for J/ip [14] and Drell-Yan 
[11] production. 

The parton densities in the gold nucleus, Fj/x u {%-,Q 2 )b,z), can be factorized into x 
and Q 2 independent nuclear density distributions, position and nuclear-number independent 



nucleoli parton densities, and a shadowing function, S^, S (A, x, Q 2 , b, z), that describes the 
modification of the nuclear parton distributions in position and momentum space. The first 
subscript on the shadowing function, P, refers to the shadowing parameterization while the 
second, S, to the spatial dependence. Most available shadowing parameterizations ignore 
effects in deuterium so that Fj/ d depends only on the deuterium density distribution and the 
nucleon parton densities. We account for the proton and neutron numbers of both nuclei. 
Then [14] 

Fi/ d (x, Q 2 , s, z) = p d (s, z)f i/N (x, Q 2 ) (4) 
F j/Au (x,Q 2 , \b- s\,z') = p Au (\b- s\,z')S 3 P:S (Au,x,Q 2 , \b - s|, z')f j/N (x,Q 2 ) (5) 

where fi/N(%,Q 2 ) is the nucleon parton density. In the absence of nuclear modifications, 
5p S = 1. The nucleon density distribution of the gold nucleus is assumed to be a Woods- 
Saxon with i?Au = 6.38 fm [15]. We use the Hulthen wave function [16] to calculate the 
deuteron density distribution. The densities are normalized so that / cPsdzpA(s, z) = A. 
We use the MRST LO parton distributions [17] for isolated nucleons and take Q 2 = p?. 

We have chosen two parameterizations of nuclear shadowing which cover extremes of 
gluon shadowing at low x. The Eskola et al. parameterization, EKS98, is based on the 
GRV LO [18] parton densities. At the minimum scale, Q%, valence quark shadowing is 
identical for u and d quarks. Likewise, u, d and s shadowing is identical at Q\. Even 
though the light quark shadowing ratios are not constrained to be equal at higher scales, 
the differences between them are small. Shadowing of the heavier flavor sea, s and higher, 
is calculated separately at Q\. The shadowing ratios for each parton type are evolved to LO 
for 2.25 < Q 2 < 10 4 GeV and are valid for x > 1CT 6 [19,20]. Interpolation in nuclear mass 
number allows results to be obtained for any input A. The parameterizations by Frankfurt, 
Guzey and Strikman combine Gribov theory with hard diffraction [21]. They are based on 
the CTEQ5M [22] parton densities and evolve each parton species separately to NLO for 
4 < Q 2 < 10 4 GeV. Although the x range is 10~ 5 < x < 0.95, the sea quark and gluon ratios 
are unity for x > 0.2. The EKS98 valence quark shadowing ratios are used as input since 
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Gribov theory does not predict valence shadowing. The parameterizations are available for 
four different values of A: 16, 40, 110 and 206. We use A = 206 for the gold nucleus with 
the parameterization that gives the strongest gluon shadowing for the largest contrast to 
EKS98, denoted FGS1 here. 

We now turn to the spatial dependence of the shadowing. Since some qualitative spatial 
dependence has been observed [23] but the exact behavior is unknown, we have tried two 
different parameterizations for inhomogeneous shadowing in d+Au collisions [11,24-26]. The 
first, S 3 P , ws , assumes that shadowing is proportional to the local density, pa{t), 



S^ ws (A,x,Q 2 ,s,z) = l + N ws [Sl(A,x,Q 2 ) - 1]^ , (6) 

Pa{0) 

where r = y/s 2 + z 2 , Pa(0) is the central density and iV ws is chosen so that 
(1/A) J d 2 sdzpA(s, z)<Sp ws = Sp- When r Ra, the nucleons behave as free particles 
while, at the center of the nucleus, the modifications are larger than S P . 

If, instead, shadowing stems from multiple interactions of the incident parton [27], 
parton-parton interactions are spread longitudinally over the coherence length, l c = 1/2itinx, 
where is the nucleon mass [28]. For x < 0.016, l c > Ra for any A and the incident parton 
interacts coherently with all the target partons in its path so that 

Sl !p (A, x, Q 2 , s, z) = l + N p [Sl(A, x, Q 2 ) - " ^ 

J azpA{v, z) 

The integral over z includes the material traversed by the incident nucleon. The normal- 
ization requires (1/A) J cPscIzpa(s, z)S P = Sp with N p > N ws . At large x, l c Ra and 
shadowing is proportional to the local density, Eq. (6). 

While there are three homogeneous FGS parameterizations, only two inhomogeneous 
parameterizations are provided. No spatial dependence is given for FGS1, the case with 
the strongest gluon shadowing. We have checked the available dependencies against those 
calculated using Sp GS1)WS and S FGSl p and found that, at similar values of the homogeneous 
shadowing ratios, Sp GS1 is quite compatible with the available FGS inhomogeneous param- 
eterizations. Therefore, to characterize the spatial dependence of FGS1, we use S FGSl p . 
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The fragmentation functions, D h / k (z c ,Q 2 ), describe the production of hadron h from 
parton k with z c = Ph/Pk- The produced partons are fragmented into charged pions, kaons 
and protons using the LO KKP fragmentation functions [29] , fit to e + e~ data. The final-state 
hadrons are assumed to be produced pairwise so that n = (n + + 7r~)/2, K = (K + + K~)/2, 
and p = (p + p)/2. The equality of p and p production obviously does not describe low 
energy hadropro duct ion well. At higher energies, however, the approximation that p = p 
may be more reasonable. The produced hadrons follow the parent parton direction. The 
minimum Q 2 in the KKP fragmentation functions is Qp ro = 2 GeV 2 , similar to but somewhat 
lower than the minimum Q 2 of the shadowing parameterizations. Thus the minimum pt 
of our calculations is \/2 GeV. We assume the same scale in the parton densities and the 
fragmentation functions, Q 2 = Qp r = p 2 T . A larger scale, p^/z 2 , is sometimes used in the 
parton densities but where z c is large, as is the case here, changing the scale does not 
significantly alter the calculated ratios. 

The largest contribution to the total final-state charged particle production is from the 
charged pions, followed by the kaons. The proton contribution is the smallest even though, 
in d+Au collisions at RHIC, (p + p)/h^ 0.24 ± 0.02 where h = h + + br for 2 < p T < 3 
GeV, independent of centrality [30]. The d+Au result is similar to that from pp, 0.21 ±0.01 
[30] . The discrepancy between the RHIC d+Au and pp results and the extrapolation from 
e + e _ is due to the poor knowledge of the fragmentation functions at large z c . 

We have calculated the pt distributions for final-state charged pions, kaons and pro- 
tons/antiprotons separately as well as the sum of all charged particles. For each final-state 
hadron, we determine the fraction of the total from produced quarks, antiquarks and gluons. 
In the central rj bin for pion production, gluons are produced almost equally in the gg — > gg 
and qg — > qg channels. The qg channel is somewhat larger for > 5 GeV. There is a 
negligible contribution from qq — > gg. Pion production by quarks and antiquarks proceeds 
mainly through the qg — > qg channel for quarks and qg — > qg for antiquarks. The next 
largest contribution to pion production by quarks are the qq' — > qq' and qq — * qq channels 
which are of very similar strength, followed by qq — > qq and gg — > qq with a negligible 
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contribution from gg — > q'q'. The contributions to antiquark production after qg — > qg are 
gg' — > qq', qq — > qq and gg — > gg, followed by smaller contributions from gg — > qq and 
gg — > q'q'. Similar results are found for kaon and proton production. However, the proton 
distributions fall off more steeply with px- 

The relative contributions from the production channels remain similar as rapidity in- 
creases. The most important change is in gluon production where the qg — > qg channel grows 
more dominant, finally becoming larger than the gg — > gg channel for all px for 77 = 3.2. 
Indeed, at the most forward rapidity, the gg — > gg channel becomes comparable to the qg 
channel at pt ~ 8 GeV. This may seem counterintuitive since the ion x<i value decreases as 77 
grows, increasing the gluon density. However, the deuteron x 1 value increases more rapidly 
and, at large 77, we are in a region where the deuteron gluon density is dropping steeply 
while the quark density, particularly that of the valence quarks, is still significant. Thus the 
qg channel is more important than the gg channel at large 77, particularly when px and X\ 
are large. Also, at high 77 and px, Pt > 7.5 GeV, antiquarks are predominantly produced by 
valence quark induced processes since these are large at high x. 

Because we begin to approach the edge of phase space with increasing 77, the px distri- 
butions steepen, especially for antiquark and gluon production. Quark production, which 
includes the valence contribution, dominant at high px and 77, remains harder overall. Thus 
quark production will come to dominate all final-state hadron production. This effect, in- 
creasingly important at high p T and 77, is reflected in the relative contributions to pion, kaon 
and proton production by quarks, antiquarks and gluons. At < 0.2, pion production 
is dominated by produced gluons up to pr > 9 GeV where pion production by quarks be- 
comes larger. Gluon production of kaons is rather small, similar to the quark contribution 
at pt ~ Qfi-o Du t dropping below the antiquark contribution at px ~ 3.5 GeV. Quark pro- 
duction is most important for protons at px > 3.5 GeV. As r\ increases, quark production of 
final-state hadrons becomes increasingly dominant. Already at 77 = 1, more than half of all 
kaons and protons are produced by quarks for px > 2.5 GeV. At higher 77, antiquarks and 
gluons make negligible contributions to low px kaon and proton production at 77 = 2.2 and 
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3.2. Quarks also dominate pion production for pr > 6.5, 3 and 1.5 GeV with r] — 1, 2.2 and 
3.2 respectively. 

We have calculated the average Au ion momentum fraction, (x 2 ), and the average 
deuteron momentum fraction, (xi), for \[2 < p T < 12 GeV. However, the largest acces- 
sible pt decreases to 9.5 GeV at rj = 3.2 due to phase space. The results are shown in 
Table I. Since these are average x 2 values, the actual x 2 for each event can be smaller 
or larger than these averages. The minimum and maximum (x 2 ) correspond to the lowest 
and highest pt values respectively. Both the minimum and maximum values decrease as 
r) increases so that the minimum (x 2 ) is reached at i] = 3.2. However, the maximum (x 2 ) 
increases relative to more central i] values due to the reduction of phase space at high pr- 
Note that as 9 — > 0, (x 2 ) — > I. The averages are not very sensitive to changes in the parton 
densities or the choice of factorization, renormalization or fragmentation scales. 

The total hadron yield closely follows that of the pions. There is little variation of 
(x 2 ) between hadron species although the proton averages are generally somewhat smaller 
than those of the mesons. A small difference between the partonic contributions to hadron 
production can be attributed to the behavior of the parton distribution functions in the 
various production channels. A set of LO parton densities derived including the GLRMQ 
recombination terms at low x found deviations from normal DGLAP evolution at x < 10 -3 
for the proton [31]. Thus one may question whether saturation effects can be at work here 
when (x 2 ) = 0.035 at rj = 3.2. 

While the average Au momentum fraction is decreasing with centrality, the average 
deuteron momentum fraction, x±, is increasing. Note that the maximum (xi) at rj = 3.2 
approaches 1, indicative of the edge of phase space. The average z c in the fragmentation 
functions is large, ps 0.5 at midrapidity, and increasing with r] and p^- The fragmentation 
functions are best determined for smaller z c so that the high z c fragmentation functions 
are unreliable, especially for baryon production. Modeling of high pt and high 77 hadron 
production thus contains large theoretical uncertainties due to the fragmentation functions. 
There is more variation in (z c ) due to parton type than in (x 2 ). The proton (z c ) tends to 
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be somewhat smaller than for pions or kaons. The tabulated values are for total charged 
hadrons. 

We now compare the ratios, i?dAu, calculated for the two homogeneous shadowing pa- 
rameterizations, to the BRAHMS data [8]. The EKS98 results are shown for each 77 interval 
in Fig. 1 while those employing the FGS1 parameterization are shown in Fig. 2. We show 
the results for charged pions (dashed), charged kaons (dot-dashed) and protons/antiprotons 
(dotted) separately. The solid curves give the total charged hadron result. At midrapid- 
ity, where (#2) is relatively large, the two parameterizations give rather similar results. As 
pointed out in Ref. [32], the difference between the kaon and proton ratios is due to isospin 
effects. As long as pion production is dominated by gluons, it is essentially independent of 
isospin. The ratio is greater than unity but smaller than the BRAHMS result at midrapidity. 
Including p T broadening would increase the \rj\ < 0.2 ratio. 

At 77 = 1 and low p?, the ratio is less than unity for both parameterizations but the 
stronger gluon shadowing in the FGS1 parameterization reduces -RdAu to ~ 0.8 for p T = \[2 
GeV relative to ~ 0.9 for EKS98. At pr ~ 2.5 GeV, i?dAu rises above unity again. At 
higher rapidities, i?dAu decreases at low pr but does not rise as far above unity at higher pr 
until, at rj = 3.2, the total charged hadron ratio is less than unity for all p T . The EKS98 
parameterization tends to underestimate the data for all but the most central rapidities, 
see Fig. 1. The FGS1 parameterization, on the other hand, agrees rather well with the 
central 77 data, Fig. 2(a), and lies within the errors of the most peripheral bins for p T > \/2 
GeV, Fig. 2(c) and (d). However, the total charged hadron data at 77 = 1 are somewhat 
underestimated by the FGS1 parameterization, Fig. 2(b). 

In the most central bins, the ratio for the total charged hadrons closely follows i?dAu 
for the pions. At higher 77, the kaon contribution becomes more important, causing the 
total to be closer to the average of the pion and kaon results. The proton contribution, 
on the other hand, remains small, even at i] = 3.2, while one may expect that, in reality, 
proton production would be more important at large rapidity as the fragmentation region 
is approached. However, this effect cannot be accounted for by standard parameterizations 
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of the fragmentation functions. Since BRAHMS measures the negative charged hadron 
distribution, h~ , at r\ = 2.2 and 3.2, only the antiprotons can contribute to -RdAu- 

We point out that because x is not very small for px larger than a few GeV and the 
shadowing ratios are also not far from unity in this region, -RdAu, the ratios of d+Au relative 
to pp, are driven by isospin rather than shadowing. This is obvious from the very similar 
behavior of the EKS98 and FGS1 ratios seen in Figs. 1 and 2 for pt > 5 GeV. 

Figure 3 illustrates the centrality dependence using S FGSlp . We compare the central-to- 
peripheral ratio, Rcp, which should be less sensitive to isospin than -RdAu, to the BRAHMS 
data. The solid curves show the ratio of the central, (0 — 20)%, to peripheral, (60 — 80)%, bins 
for each rj region while the dashed curves show the semi-central, (30 — 50)%, to peripheral 
ratios. Our calculations assume exact impact parameter cuts while, experimentally, impact 
parameter is poorly measured on an event-by-event basis in d+Au collisions. We note that 
both of the inhomogeneous EKS98 results are much weaker than those in Fig. 3 and are not 
shown here. 

In central collisions, with small impact parameter, inhomogeneous shadowing is stronger 
than the homogeneous result. The larger the homogeneous shadowing effect, the larger the 
difference between Sp and Sp s . Thus R G p is a stronger function of impact parameter for the 
FGS1 parameterization since it has larger homogeneous shadowing at small x. The ratios 
with Seks g underestimate the centrality dependence considerably and are not shown. Note 
that i?cp approaches unity at large px since the difference between Sp and Sp s decreases 
as x increases and the shadowing effect becomes small. 

Since S FGSlp (b) is a rather smooth function of impact parameter, the b dependence of 
Rcp is not very strong. The fluctuations in R C p for S FGS1 , especially notable at central 
rapidity, are due to the discrete steps of Tj^ u (r) in the integration over the spatial coordinates. 
These fluctuations are absent for Sp ws since PAu( r ) is a smooth function. 

The agreement with the data is reasonable at central rj, see Fig. 3(a) and (b). The 
trends of the impact parameter dependence are similar to the data at low px- The semi- 
central-to-peripheral ratio is similar to the central-to-peripheral ratio in the most central 
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rapidity bin while in the more peripheral bins, the central-to-peripheral ratio has a stronger 
Pt dependence. However, the increase in (x 2 ) with p T results in the strong growth of Rqp 
with pt at forward 77. The resulting curvature of the calculated ratio is faster than the data. 
The magnitude of Rqp at low pt is also underestimated. Since the position dependence of 
inhomogeneous shadowing is not well understood, the poorer agreement with the centrality- 
dependent data in Fig. 3 compared to the minimum bias results in Figs. 1 and 2 is not 
surprising. These data could be used to tune the position dependence of shadowing. 

In summary, we find that the suppression factor, -RdAu, calculated with leading- twist 
shadowing, especially employing the FGS1 parameterization, agrees moderately well with 
the BRAHMS data. These calculations imply that saturation effects may not play a dom- 
inant role in the forward region at RHIC, as suggested in other recent work [33,34]. Our 
calculations of Rqp show a stronger p T dependence than that suggested by BRAHMS, likely 
due to insufficient data on the impact parameter dependence of nuclear shadowing. 

We thank K.J. Eskola and V. Guzey for providing the shadowing routines and J. Gonza- 
lez, S. R. Klein, M. Murray, M. Strikman and W. Vogelsang for discussions. This work was 
supported in part by the Division of Nuclear Physics of the Office of High Energy and Nuclear 
Physics of the U. S. Department of Energy under Contract Number DE-AC03-76SF0098. 
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TABLES 



(v) 


<x 2 > 

min 


max 


min 


max 


(Zc) 

min 


max 


« 


0.07 


0.22 


0.07 


0.22 


0.52 


0.64 


1 


0.055 


0.18 


0.1 


0.33 


0.54 


0.68 


2.2 


0.042 


0.14 


0.17 


0.62 


0.55 


0.81 


3.2 


0.035 


0.23 


0.32 


0.95 


0.64 


0.96 



TABLE I. The average values of the Au and d momentum fractions, {x2) and (x\) respectively, 
as well as the average fraction of the final-state parton momentum transferred to the hadron, (z c ), 
in the four BRAHMS pseudorapidity intervals. The minimum values correspond to px ~ \/2 GeV 
while the maximum corresponds to px = 12 GeV for the first three r\ bins and 10 GeV for the most 
forward rj bin. 
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FIG. 1. -RdAu f° r charged pions (dashed) and kaons (dot-dashed) as well as protons and an- 
tiprotons (dotted) and the sum over all charged hadrons (solid) for deuteron-gold collisions at 
y / Sjv r /v = 200 GeV as a function of pr- The results for homogeneous shadowing with the EKS98 
parameterization are compared to the minimum bias BRAHMS data [8] in the following r\ bins: 
(a) \r)\ < 0.2; (b) r? = 1; (c) r? = 2.2 and (d) rj = 3.2. 
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FIG. 2. i?dAu for charged pions (dashed) and kaons (dot-dashed) as well as protons and an- 
tiprotons (dotted) and the sum over all charged hadrons (solid) for deuteron-gold collisions at 
V$nn = 200 GeV as a function of pr- The results for homogeneous shadowing with the FGS1 
parameterization are compared to the minimum bias BRAHMS data [8] in the following 77 bins: 
(a) \q\ < 0.2; (b) r? = 1; (c) r? = 2.2 and (d) r? = 3.2. 
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FIG. 3. i?cp for charged hadrons in deuteron-gold collisions at \/Snn = 200 GeV as a function 
of pt- The results for Sfgsi.p are compared to the BRAHMS data [8] in the following rj bins: (a) 
\v\ — 0-2; (b) r7 = l; (c) rj = 2.2 and (d) rj = 3.2. The calculated ratios of the most central and 
semi-central to peripheral collisions are shown in the solid and dashed curves, respectively. The 
BRAHMS data are given by the open circles (most central) and diamonds (semi-central). 
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